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INTRODUCTION 


The purpose of this paper is to examine and correlate the 
various physical and chemical relations proposed to explain the phenom- 
ena of detonation, in particular, the detonation of solid explosives. 

In order to assure a reasonably comprehensive treatise, it is desirable 
at the outset to devote some attention to the correction of certain 
widespread misconceptions regarding explosives and their behavior. Two 
distinct types of explosion processes, deflagration and detonation, are 
evidenced by materials referred to as explosives. It has been common 
practice, especially in the early literature on the subject, to distin- 
guish between the two on the basis of chemical composition of the 
explosive and its intended application. Deflagration is invariably 
associated with "low" explosives or propellants and the continual employ- 
ment of certain explosives such as nitrocellulose and black powder for 
their propulsive effect has led to the erroneous association of chemical 
compositions with the type of behavior. Similarly the detonation process 
is associated with primaen and "high" explosives and likewise is con- 
sidered to be the sole type of behavior evidenced by other chemical 
substances. Unfortunately the above classification is not borne out by 
facts. On the contrary, experimental results indicate that all explosives 
are capable of exhibiting both modes of behavior. In view of this it is 
necessary to establish some criterion other than chemical composition to 
enable one to distinguish clearly between the two processes. 

The deflagration process, in its broadest sense, can be con- 
sidered a rapid burning in which the essential rate-determining factor 


is the speed of the chemical reaction, That is, the rate at which the 
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deflagration proceeds is related directly to the kinetics of the reaction 
taking place in the burning zone. 

The propagation of deflagration is closely associated with 
transport mechanisms; either the transport of heat or of atoms and free 
radicals acting as chain-carriers or a combination of the two. As long 
ag mechanisms such as these control the propagation, there is established 
a limiting velocity beyond which the rate of deflagration cannot pass. 
Extensive investigations have been carried out in connection with the 
process of deflagration in a gaseous medium and several excellent treat- 
ments are available (1,2). 

In the case of detonation, the characteristic feature distin- 
guishing the process is the appearance of a pressure wave propagating at 
a high velocity. The rate at which detonation proceeds is calculable by 
means of hydrodynamic and thermodynamic equations without the necessity 
of reaction-kinetic data. The application of the classical hydrodynamic- 
thermodynamic theory requires no assumptions regarding chemical mechanisms 
other than the postulate that the reactions occur at such a rate so as 
not to be a limiting factor. 

This theory, which is reviewed in detail, is based upon 
expressions for the continuity of mass, momentum and energy across a 
so-called detonation front and upon the assumption that the stable 
detonation velocity is the minimum one compatible with the conservation 
equations. When supplemented by an equation of state and thermal data 
for the products of reaction, the theory permits the calculation of 
detonation velocity, temperature and pressure, as well as other parameters 


of interest. 
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INITIATION AND PROPAGATION OF DETONATION 





The currently accepted hydrodynamic detonation theory, which 
will be discussed in a later section, unfortunately provides no informa- 
tion concerning the mechanism by which a stable detonation wave arises. 
It is well established by experience that detonation can be brought about 
by a wide variety of conditions, including mechanical shock, high 
temperature, friction, and adiabatic compression of interstitial gas. 
Only recently, however, have successful attempts been made towards the 
establishment of what might be considered a generalized theory of ini- 
tiation. The work of S. J. Jacobs and D. P. MacDougall (3) at the 
Explosives Research Laboratory of the N.D.R.C. involving the controlled 
deflagration of readily detonable substances and that of F. P. Bowden 
and his associates (4) on the isolation of factors responsible for 
thermal ignition provided a basis upon which such a theory could be 
proposed. The approach adopted by both G. B. Kistiakowsky (5) and A. R. 
Ubbelohde (5) is that of a preliminary deflagration evolving discontinu- 
ously into a detonation. 

Consider a mass of explosive particles in which a small 
peripherial deflagration has been started. As the deflagration spreads 
inward, pressure gradients develop due to the restricted flow of the 
gaseous products through the interstices between the particles. As the 
interior pressure rises, the speed of the deflagration increases, result- 
ing in a faster decomposition and a further increase in pressure. This 
constantly increasing deflagration may evolve into a stationary combus- 
tion in which the flow of heat and material away from the reaction zone 


compensates for that produced by the reactions, or the deflagration rate 
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may continue to increase until all the material is consumed provided the 
mass of explosive involved is less than the "critical size". The concept 
of "critical size" for a particular explosive is a function of many 
variables including shape (volume-area ratio), particle size, crystal 
form for polymorphic substances, and density. It is probably best 
defined, although somewhat nebulously, as a limiting criterion, that is, 
the critical size is that size at which the rate of decomposition under- 
goes a discontinuity. The discontinuity may result in a partial or 
complete extinguishing of deflagration, with perhaps a re-initiation at 
an adjacent spot in the substance or it may be evidenced by the phenomenon 
of detonation. To see how the non-stationary deflagration might evolve 
into a detonation we must consider the flow of gases from a center of 
deflagration in more detail. As each crystal decomposes it acts as a 
center from which the gas molecules stream outwards, thus forming a 
pressure wavelet. If a number of these wavelets can integrate to give 
a streaming velocity corresponding to that accompanying a stable 
detonation wave for the system, we have conditions identical to those 
observed in a detonation. The type of flow just postulated has a very 
close analogy in the field of gaseous combustion. A column of gas acted 
upon by a piston subject to acceleration is frequently cited as an example 
of the mechanism of shock wave formation in gases. Further discussion of 
this mechanism is included in the following section. 

Assuming a shock wave has developed it is of interest to con- 
sider how its passage through the explosive can initiate the decomposi- 
tion. The details of the mechanism by which a shock with its discontin- 


uous rise in pressure and non-isentropic temperature effects activate 
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the explosive molecules is still the subject of considerable specula- 
tion. 

It has been proposed by Ubbelohde (5) that the most probable 
mechanism for the propagation of detonation is the collision of high 
speed atoms, radicals or molecules from behind the wave front with the 
quiescent layer of molecules just ahead of the front. Particles behind 
the detonation front are assumed to have an unusually high activation 
efficiency due to high random velocities at the detonation temperature 
and a high directed velocity resulting from the hydrodynamic conditions 
of propagation of a shock wave. Impact of a high speed particle is 
visualized as disrupting an unreacted molecule before the kinetic energy 
has had time to be distributed among all the vibrational degrees of 
freedom. If dissociation occurs after the redistribution of energy, 
then such a dissociation must be attributed to high temperature in the 
ordinary sense. Since redistribution requires a small or zero activation 
energy, whereas dissociation requires an activation energy about equal 
to the strength of the bond broken, it is probable that redistribution 
of the kinetic energy will occur before dissociation in most cases. 

It has been suggested by Wendlandt (6), Eyring (7) and others 
that the mechanism of chemical activation depends primarily on the high 
temperature produced by the shock wave. Although considerable 
uncertainty exists in regard to reactions at the extreme temperatures 
and pressures of detonation, it has been shown (8) that the temperature 
coefficient of the energy-releasing process is much smaller than the 
usual Arrhenius value. This, plus the fact that it is unlikely that 


sufficient time exists in the activation process for the vibrational 
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and rotational quanta to attain the full Maxwell-Boltzman equilibriun, 
indicates that the temperature effect alone probably does not completely 
account for the activation. 

Still another theory has been proposed suggesting that the 
pressure in the wave front causes a bending of the molecules placing 
reactive groups in close proximity to each other, whereupon the 
decomposition occurs. It is felt that this theory is the least satis-—- 
factory of those presented since Bridgeman (9) has shown that high 
hydrostatic pressure alone has no effect on explosives. This theory is 
also untenable from the standpoint of chemical mechanism since practically 
every conceivable first step in the decomposition of organic nitro 
explosives is endothermic and the pressure activation concept requires 
an exothermic initial step. 

Thus at the present time it appears most reasonable to 
attribute the propagation of detonation to the energy possessed by the 
shock wave as a whole rather than to attempt a more detailed and perhaps 


erroneous description of the phenomenon. 
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SHOCK WAVES 


The theoretical basis for the development of a theory of 
detonation stems from theory of shock waves proposed by Hugoniot (10) 
in the nineteenth century. In mathematical language he formulated the 
difference between shock waves, or compressional waves of finite 
amplitude, and weak acoustic waves. The former were found to be 
propagated with velocities which are greater than the acoustic velocity 
of the medium, the more so the higher the peak pressure in the crest of 
the wave. The medium itself is moving in the shocks, the wave velocity 
being slightly less than this mass velocity plus the acoustic velocity. 
Regardless of the original shape of the (finite) compressional wave in 
the medium, a shock front is formed eventually as the wave propagates 
through the medium. This means that pressure, density, and temperature 
suffer discontinuous increases in the front of the wave. In the course 
of its propagation through a medium a pure shock wave decays, this 
process being superimposed upon the normal weakening of waves due to the 
geometric factor of divergence, unless mechanical energy is constantly 
fed into the shock wave. In other words, the mechanical energy of shock 
waves is being degraded and is then found as the heat energy of the 
medium after the passage of the wave. 

Since a detonation wave is basically a self-sustaining shock 
wave, a brief qualitative discussion of shock wave formation appears 
desirable. It is well established that in a compressional wave of 
finite amplitude the velocity of a portion of a compressional pulse is 
directly proportional to its amplitude. Those portions of the pulse 


with the greatest amplitude will travel faster than portions of lower 
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amplitude, theretore changing the shape of the pulse as indicated in 


Figure l. 


b b b 
; , a ws ta 
Figure 1 
Formation of a Shock Front Pp >P, 


Becker's treatment (11) of the origin of a shock wave in a 
gas is also of interest. Let a piston in a long tube be given a sudden 
small velocity. The gas in front of the piston will be compressed and 
a pressure wave will move ahead of the piston with the velocity of sound 
in the medium. The gas in front of the piston and behind the pressure 
wave will be moving with the velocity of the piston. If the velocity 
of the piston is increased by another small increment, a new pressure 
wave arises that moves somewhat faster than the first, since it moves in 
the gas already heated by the forezoing compression and moving with the 
velocity of the piston. kach succeeding wave will move more rapidly 
than the preceding one and will ultimately overtake the initial one to 
form a composite wave with an infinitely steep front; i.e., a discon- 
tinuity in pressure and temperature will be formed. 

The qualitative argument can be made more explicit by the 
Riemann formulation of the equations of hydrodynamics. Such a treatment 
is presented by Lamb (12) and the results confirm the conclusions reached 
above. 


By reversing the argument given for compression waves, it is 
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apparent tnat a rarefaction wave of finite amplitude, in which the later 
portions of the wave are regions of lower pressure, will spread out as 
it advances and a discontinuity will not arise. 

Returning to the shock front, which for the purposes of this 
discussion can be treated as a mathematical plane, it becomes desirable 
to establish equations relating the conditions on both sides of the 
discontinuity. The equations for the conservation of mass, momentum and 
total energy have been employed by Rankine (13) and Hugoniot (10) in the 
formulation of three conditions relating the pressure, specific volume 
and particle velocity benind a discontinuity moving with a velocity D to 
those properties in the undisturbed medium ahead of the discontinuity. 

Consider a shock front moving into a region in which the 


particle velocity is u,, the density Po and the pressure P, from a 
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region in which the same properties are denoted by u,,( and P. 





shock front 
Figure 2 
Conditions at a Shock Front Moving with Velocity D. 


In a time dt a mass of fluid C_(D-ug)dt will enter unit cross-section 
area of the front and a mass C1 (D-u, )dt will leave the front (see 

Figure 2). As dt is made very small so that the layers on either side 
are infinitesimal, the mass flow into the front must equal the mass flow 


away from it. Thus an expression for the conservation of mass may be 


written: 
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C,(D-u,) = ,(D-u) (1) 


The change in momentum associated with the mass flow across 
the front must equal the net force per unit area. Thus for the con- 


servation of menentum: 


Po (D-ug) (uz -u,) = P)-P, (2) 


Finally the conservation of energy requires that the net work 
doneby the pressures P) and P, equals the increase in kinetic and 
potential energy. The work done by each pressure is Pudt for unit area, 
the kinetic energy is 4 u* and the internal energy is E per unit mass. 


Therefore, 


PiUj)-Poy = (%(D-Ug) Hy -Eo+4(uf-ué) (3) 


The above equations may be rewritten in the form below by 


substitution of V = 1/6 and some algebraic manipulation. 


D= ugt¥, J (Py-P,)/W-V1 (4) 
Uy = uot " (Py -P,)(V_-V}) (5) 
E)-E, =a) (PSAP, )(V6-¥,) (6) 


The three equations above are often referred to as the Rankine- 
Hugoniot conditions. 

Before proceeding to the consideration of detonation processes, 
it is desirable to consider Equation 6 in more detail and to investigate 
the thermodynamic properties of the shock wave. Knowledge of the 
equation of state and the heat capacity enables one to calculate 


E)-E, as a function’of P, and V, so that Equation 6 becomes the Rankine- 
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Hugoniot equation relating the pressure and volume in a compressional 

wave of finite amplitude for given initial conditions. The relation 
describes the change in internal energy before and after compression by 
the shock wave as a function of the pressures and specific volumes. 

This change is adiabatic, since no heat is exchanged with the surroundings 
but it is not isentropic. The non-isentropic nature of a shock (or 
detonation) wave was not appreciated by many of the early investigators 
and led to a conflict with the energy principle. A method of determin- 
ing the increase in entropy occurring across a shock front is discussed 


along with other detonation zone parameters in a later section. 
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DETONATION WAVES 


The existence of detonation waves was originally reported by 
Berthelot and Vieille (14) in connection with studies of gaseous 
explosions. Such a wave differs from a shock wave in two respects. In 
the first place, the wave is sustained and propagated through the 
mechanism of a chemical reaction supplying energy. Therefore, the 
chemical energy released must be included in the calculations of the 
Rankine-Hugoniot conditions. The second difference is that the velocity 
of the wave is determined by the thermodynamic state of the material 
immediately behind the wave front rather than by the motion of an 
external boundary surface as is the case in pure shock waves. The three 
conditions relating properties of the material on either side of the 
detonation front in terms of the conservation of mass, momentum, and 
energy are no longer sufficient to determine the unknown quantities 
behind the front in addition to the velocity of the front. The problem 
differs from that considered in the previous section in that there is 
no controllable external condition, such as the velocity of a piston, 
available to determine one of the variables behind the front. 
The Chapman-Jouguet Condition 

An additional relation is required to make possible a unique 
solution of the hydrodynamic equations. Chapman (15) postulated a 
relation between the velocity of the detonation wave and the sonic 
velocity in the gas behind the wave. He also assumed that the mininun 
detonation velocity consistent with the Hankine-Hugoniot conditions pre- 
vailed a a steady detonation. Although his final results were in error 


due to assuming that isentropic conditions prevailed at a discontinuity, 
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the basic concept was ide:itical to the mathematical expression published 
by Jouguet (16) two years later. The condition now bears the name 
Chapman-Jouguet and postulates 


D = C) +; (7) 


where uj, is the material velocity and c) the sonic velocity in the 
reaction products behind the wave front. 

Although a valid theoretical proof of the necessity of this 
condition is not yet available, there is little doubt that it correctly 
describes the observed phenomenon. Several plausible arguments have 
been advanced which aid in the acceptance of this condition. The general 
form of the pressure-volume curve representing possible states of the 
products after detonation is shown in Figure 3. This curve is specified 
by the Rankine=Hugoniot conditions, the thermochemical data concerning 


the explosive and the equation of state. 





Figure 3 
Generalized Rankine-Hugoniot Curve 


The point A represents the initial conditions existing in the 


explosive before the passage of the detonation wave. If the initial 








=a 


material velocity is zero, as is invariably the case with solid 
explosives, Equation 4 gives the detonation velocity in terms of the 
initial and final states. 


D=Vo ff (P1-Po)/Vo-V1 (8) 


From Figure 3 it is evident that this relation can be written as 


D=V,y tan 6 (9) 


where @ is the angle between the system of lines passing through point A 
and the (negative) V axis. For any value of 9 greater than 0, there 
will be two distinct final states (C and G) corresponding to the same 
value of D. For 6 = 0, the line AB is tangent to the curve at point B 
and is the final state corresponding to the minimum value of D. 

If the final state of the products corresponds to a point C 
on the R-H curve higher than the point B, it can be shown (see Appendix A) 
that the velocity of the detonation wave is less than u+c and a rare- 
faction wave travelling at sonic velocity with respect to the products 
would move faster than the detonation and hence engulf and destroy it. 
(A rarefaction wave will arise from any process that abstracts energy 
from the reaction products such as heat loss or turbulence. ) 

At the point B, the velocity of the detonation wave is equal 
to the velocity of sound in the reaction products plus the mass velocity 
of those products, thus a rarefaction wave proceeding with its maximum 
velocity will be unable to overtake the detonation wave. The mathe- 
matical proof that D = ceu at point B follows that given by Kistiakowsky 
and Wilson (17). From Equations 4 and 5 the detonation velocity is given 


by 
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V 
D = 2 uy (10) 
ver 1. 





The velocity of sound in a gas is given by the Ae esis (18) 
OP 
caf (aP7 ape = v(-3el (11) 


If there is a point on the Rankine-Hugoniot (R-H) curve 
having the same slope as the corresponding point on the isentropic curve 


(PV curve with constant entropy), then at that point 


QF) «2 /Zt 








where the subscript r refers to the Rankine~Hugoniot curve. Differ- 
entiation of the Rankine-Hugoniot relation, Equation 6, holding P, and 


Vo constant yields 


dE, = - 2(Po+P))dVy+3(V-V, )dPy (12) 


Now, from the first law of thermodynamics, dE = TdS-PdV, thus Equation 


12 can be written 
dS} dP) Je al 
Ty ioee= | F BV oMV eae) ~ a 1 
1 2). 5 ( ° 1) fan 7 ard (13) 
It follows that the slope of the R=-H curve through a point at which 


dS, = 0 can be expressed as shown in Equation 14 below 
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This is also the slope of the line AB which is tangent to the R-H curve 
at point B. Now the detonation velocity at any point on the R-H curve 


is given by Equation 4 and is 


-p.\3 
D=V, Po2 (4) 
Vor¥) 
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Therefore, at point B 
ap, \2 
D=Vy [-et} = 22 ¢, (by Eq. 1) (25) 
Combining this with Equation 10 yields 


The argument which is used to exclude points on the R-H curve 
below B is based on the entropy of the products. Consider the two points 
C and G lying on the same straight line from A and corresponding to the 
same detonation velocity D. The entropy of products at C is greater than 
at G since the transition from G to C is equivalent to passage from an 
initial state G across a shock wave of velocity D to the final state C. 
Since there is a definite increase in entropy accompanying the passage 
of a shock wave the entropy at C will be greater than at G, therefore 
the upper point will be a thermodynamically more probable state than G. 
This is true for all angles of 6>6,. 

Thus points above 8B are more probable than points below but 
points above B are not stable due to the effect of the rarefaction wave. 

The question of the validity of the Chapman-Jouguet condition 
has also been discussed by von Neumann (19). By applying the conserva- 
tion Equations 4, 5 and 6 not only to the completely reacted material 
but also to each stage of the reaction, he has been able to substantiate 
the validity of the condition in certain cases. In addition to the 
Single R-H curve in Figure 3 representing the final state of the 
products, there will be a series of R-H curves corresponding to the 
intermediate stages. Since the detonation velocity is the velocity with 


which the region of chemical reaction progress, it is the same for each 
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intermediate stage. Therefore, the line from the initial state P., V, 
to each Pj, Vp) representing the condition of a given layer must have 
the same slope and must therefore be the same line for all degrees of 
completion of the reaction. Since in order to get from incomplete to 
complete reaction, the material must have traversed all intermediate 
stages, this line must touch all of the R-H curves corresponding to 
different degrees of completion. The family of R-H curves corresponding 


to the degree of completion of the chemical reaction from n = 0 to 


n= 1 has the general form shown in Figure 4. 


Reaction Complete 





Figure 4 
Family of Rankine-Hugoniot Curves fromn =O ition=l. 


The motion of a point on the diagram representing the behavior of an 
element of explosive passing through the shock front and into a finite 
reaction zone must always remain on a straight line to satisfy the 
conservation equations. In the shock wave, the representative point 
must begin on the curve n = O, at point A, must intersect the curves 
for increasing n in order and finally end somewhere on the curve n = 1. 


In order that there exist a discontinuity in pressure and temperature 
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at the shock front, there must be another intersection of the straight 
line with the curve n = 0, corresponding to point C. Thus in the shock 
front the representative point runs from A to the limiting value at the 
point C where it again meets the curve n = 0. So far no chemical 
reactions have occurred but the explosive material has been "shocked" 
into a condition where reactions can proceed. Now the representative 
point must proceed continuously from the curve n = 0 ton =1; this 
requirement immediately restricts the possible intersections to points 
onn #* ] at or above B. If the straight line intersects the curve 

n= 1 at a point above B, the tangent to the curve at that point would 
be steeper than at B, so the velocity of sound in the products would 
exceed the detonation velocity, and a rarefaction wave would engulf and 
destroy the detonation wave. Thus the intersection with the curve 

n = 1 must be at B giving a point of tangency and corresponding to Chapman- 
Jouguet condition, D = uj,+c). 

Von Neumann also discusses the case of intersecting R-H curves, 
but since this possibility does not correspond to any known physical sit- 
uation it will not be pursued further. 

By a straightforward thermodynamic procedure, Scorah (20) has 
shown that when the Chapman-Jouguet condition is assumed the ratio of 
entropy increase to available energy passing through the shock front is 
a maximum. Although this result is in accord with the fundamental 
principle of maximum degradation of energy observed in natural processes, 
it is regarded only as a further substantiation of the condition rather 


than a rigorous proof. 
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THE DETONATION ZONE 


Before proceeding with the development of the fundamental 
theory, it is advantageous to have an idea of the order of magnitude of 
the dimensions of the zones in which the process of detonation occurs. 
The essential elements of a steady one-dimensional detonation process 
can be illustrated as in Figure 5. The detonation zone is considered 


to be stationary in space and the explosive feeds into the zone from 


right to left at the detonation velocity D. 
ee Det Oo zone 






Shock Pos Vo 
Reaction Zone zone 
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Steady Plane Detonation in a Solid Explosive 


Conditions throughout the zone are considered to be steady. The material 
first passes into a plane shock zone where abrupt increases of pressure, 
temperature, mass velocity and entropy occur. This transition takes 
place within a narrow region whose length is determined by the viscosity 
and the thermal conductivity. If the dimensions of this shock zone are 
very small in comparison with other lengths under consideration, the 
shock may be treated as a mathematical plane or discontinuity without 
error. For one-dimensional motion, the Navier-Stokes equations (21) of 


motion and energy for a steady state in which the shock wave is at rest 


are: 








and 


: 2 2 Vo AT 
QE ~ 69 -4(D-u)® -Po(VorV) = k Zo 


where m is the viscosity, k is the thermal conductivity, u is the 
material velocity, T is the temperature and x is the coordinate in the 
direction perpendicular to the shock front. These differential equations 
have not as yet been solved in closed form for general values of the 
parameters. For the case of a pure shock wave (4Q ® 0) in liquid di- 
ethyl ether, Becker (22) has calculated the length of the zone to be 
6.5A° at a shock pressure of 10,000 atms. and 1.4A° at a pressure of 
100,000 atms. For a solid explosive with an assumed viscosity of 0.2 
centipose and a thermal conductivity of 104 ergs om7* sect, Eyring (23) 
calculated the length of the zone to be of molecular dimensions. Calcu- 
lations including the effects of viscosity and thermal conductivity in 
the conservation equations have shown these effects to be of importance 
only at the very front of the detonation wave, where the amount of chem- 
ical reaction is negligible. 

The distance between the points of initiation of the chemical 
reaction and the point at which the Chapman-Jouguet condition is 
applied is called the reaction zone. There are a number of methods for 
estimating the length of this zone. Since each method rests on the 
interpretation of a different phenomenon, agreement among results found 
by different methods greatly enhances the accuracy of the determination. 
The results of (1) direct observation of the luminous zone, (2) extrapo- 
lation of low-temperature rate measurements, (3) surface reaction rate 
calculations for granular explosives, and (4) studies of the building-up 


of a low velocity wave to steady state conditims, all indicate that the 
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reaction zone for INT under steady state conditions is about 4 centimeter 
in length. Although the exact figure will vary from explosive to explo- 
Sive, it is reasonably well established that the reaction zone length 
for solid explosives is between ] mm. and 1.5 cm. 

The time required for a particle to pass through the reaction 
zone can be calculated knowing the zone length, detonation velocity, and 
the material velocity. Reaction times of the order of magnitude of 


1 microsecond have been found for practically all organic explosives. 
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THE HYDRODYNAMIC THEORY 


Application of the Rankine-Hugoniot relations and the 
Chapman-Jouguet condition to an actual chemically supported detonation 
permits evaluation of the energy relation, Equation 6, in terms of heat 
of reaction and specific heat data for the reaction products. Since the 
shock zone is of molecular dimensions any contribution of heat by thermal 


conduction is neglected. The energy relation can now be written 


L 
= 1 ak 
E, 7H, = 3(P +P.) (V,-Vz) — - + rif Cydt+ 32 dV (16) 
T> a) Ty 
or introducing the mean specific heat at constant volume and assuming an 


ideal gas 


3(P,+P,)(V-Vy) = 2 + Z (T,-T,) (17) 


where Q is the heat absorbed in the reaction of M grams at the initial 
temperature T,. As indicated above the computation of conditions about 
the detonation front requires the knowledge of the following: 

1. Equation of state of the products 

2. Composition of the reaction products 

3. Heat capacity of the products 

4. Heat of reaction 
In the case of solid explosives, the selection of an equation of state 
for products at temperatures of the order of 3000 = 5000° A and at 
pressures of 50,000 to 200,000 atmospheres presents a formidable problem. 
Our present day knowledge of equilibrium constants under the above 
conditions is also quite meager. In view of this situation, it appears 


desirable to test the validity of theory under the less stringent con- 


ditions provided by gaseous explosions. 
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Gaseous Explosives 
The approach adopted in this discussion is similar to that 


applied to mixtures of hydrogen and oxygen by Lewis and Friauf (24). 
Since the equations used for ideal gases will later prove of value in 
the preliminary treatment of solid explosives, they will be fully 
developed here. The Kankine=Hugoniot conditions for conservation of 
mass and momentum, Equations 1 and 2, may be written as below, if 


Uy is assumed to be zero. 


D = De-uy (18) 
Yo Vy 

Du) . 6 < 

ic = P)-Py (19) 


Klimination of u, from the above equations yields 


2 
=p 
Py-Po * y2 (V,-Vz) (20) 


At the Chapman-Jouguet point on the R-H curve the detonation velocity 


expression Equation 4 can be written 


D = Vo ’ (Py -PS WV) = Ye r -(d4P,/aV,), (21) 


For an ideal gas, the isentropic expansion law is 


V1 
PiVy = constant 
Therefore, differentiating and substituting into Equation 21 
D = V/V} [Pi %y =~ V/V} {yay RT| /M (22) 


in which n, is the number of moles of gas per M grams of product gases. 


Since the initial pressure P, is negligible in comparison with P, in 
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Equation 20, substitution of Equation 20 into Equation 22 with the 
elimination of D yields 


_) a 
Ye 


on) 





. (23) 


This expression may now be substituted into the energy relation, 


Equation 17, again neglecting Po to give 


n,RT = 
3 — Q + C,(T,-T,) (2h) 


If the chemical reaction under consideration went quantitatively and 

the distribution of atoms in the products was independent of temperature 

and pressure, Equation 24 could be solved for the final temperature Ty 

if the variation of V1 and C, with T were known. Actually, of course, 

both ny and Q depend upon pressure and temperature. The calculation is 

therefore somewhat more tedious but the principle involved is unchanged. 
Once the final temperature is found, the detonation velocity 


can be computed from the cambination of Equations 22 and 23 


D= (841) J njRT}/ 04M (25) 


Lewis and Friauf have calculated detonation velocities for mixtures of 
hydrogen and oxygen with various excesses of both components and with 
the addition of inert gas. They allowed for the equilibria 

B Hot0, = H,0, Hp = 2H, and H,0+30. = 20H. Although the dissociation of 
QO» was not considered and no allowance was made for the excess heat 
capacity of oxygen molecules due to electronic excitation, there is no 
reason to suggest that a similar calculation undertaken today would 
differ significantly from their results. The following table shows 


their calculated values and the experimental results obtained by Dixon 
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(25) and Payman and Walls (26). 


Table I 
Detonation Velocity 
D, m/sec 

Mixture Computed Observed Deviation % 
(2H,+05 ) 2806 2819 ~O.4 
(2H5 +0 + 105 2302 2314 -0.7 
(215405) + 305 1925 1922 Ove 
eras + 505 a o2 1700 2.0 

2H2+0, ) + IN5 2378 207 -l.2 
(2iz402 + 3N5 2033 2055 -1.1 
2Hoat02) + 5No 1850 1822 Lod 
(2H2+02) + 2Ho 3354 3273 2.5 
(2Ho+02) + AHo 3627 3527 2.8 
(2Ha+02) + 6H2 3749 3532 6.1 


The agreement between the calculations and the observed measurements 
serves as an extremely good verification of the hydrodynamic-thermo- 
dynamic theory. With one exception the differences between the 
comparative values are all within the experimental errors of velocity 
measurements. The detonation velocity of a number of gaseous explosives 
of the hydrocarbon-oxygen type also have been calculated and found to 
be in good agreement with experimentally determined values. 
solid Explosives 

The extension of the fundamental theory to the detonation of 
solid explosives requires the selection of an equation of state applicable 
to gases under extreme conditions of temperature and pressure. Use of 

the ideal gas law is immediately ruled out since the experimentally 

| observed detonation velocities are quite sensitive to initial density 
and Equation 25 based on the ideal gas law shows no dependence on this 
factor. 


Another complicating factor arising with solid explosives is 








|e 


the question of the composition of the reaction products. To consider 
all of the various equilibria possible when a typical explosive contain- 
ing atoms of carbon, hydrogen, nitrogen and oxygen detonates is a 
tremendous mathematical task and cannot be justified in view of the 
uncertainty of the thermodynamic functions of the possible products. 
Justification for ignoring the more complicated equilibria involving 
free atoms, polyatomic molecules and free carbon in approximate 
computations will be presented later. 


Variation of Parameters = Abel Equation of State 





To illustrate the form assumed by the hydrodynamic equations 
and to ascertain certain qualitative aspects of the process of detona- 
tion a particularly simple equation of state will be chosen as an 
example. An equation suitable for this purpose, but of less value in 
connection with practical calculations, is the Abel equation, 

P(V- &%) = RT, with the covolume & assumed constant. 

The details of deriving the relations will for the most part 
be omitted. The derivations are quite straightforward and similar to 
those given by kyring and his associates (7). 

Evaluation of the sonic velocity appearing in the Chapman- 
Jouguet conditions, Equation 7, in terms of the equation of state yields 


(Deu, )* = V,*y - ~ (26) 
ae 





which when inserted into the Rankine=-Hugoniot relations and combinations 
of these relations gives the following five equations. (The initial 


pressure P, is negligible in comparison to P, and is omitted.) 
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Vy . vw + CYV, 


Vo 2 wT 41 (27) 
Py(Vi-ex) = &R Attyte (28) 
Cy 
Q+CyTo 23 
Ae ——— 2 
D* = (Qe0,T,) 20 eey ; (30) 
(14 &/Vo) 


a/V 
cts TEs ap (31) 


we +1 


Equations 27 to 31 furnish the detonation velocity and all other 
properties of the reaction products provided we have data concerning 
the heat of reaction, the mean heat capacity of the products, the co- 
volume and the gas constant per gram. It is apparent that neither the 
detonation velocity nor any of the properties of the products depend 
upon the rate of chemical reaction. 

| An estimate of the probable accuracy of the calculated 
properties can be made from the above equations. The final temperature 
is a function of the heat of reaction and the mean heat capacity of the 
products and in particular is independent of the covolume; thus the 
temperature should be determinable to a higher accuracy than properties 
depending on additional data. In particular the final pressure is least 


accurately obtained since, in addition to dependence upon the thermal 
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data and the covolume, it is also a function of the gas constant which 
depends upon the mean molecular weight of the products. 

The properties of the reacting explosive can be determined 
at various intermediate stages by replacing Q in Equations 27 to 31 by 
nQ where n denotes the degree of completion of the reaction. It is 
interesting to note that such a procedure is another application of the 
mechanism used by von Neumann in his attempt to justify the selection of 
the Chapman-Jouguet condition based upon the family of Ranxine-Hugoniot 
curves. The variation of the properties, temperature, pressure, ratio 
of c+u to D and entropy of a typical explosive as the reaction proceeds 
is indicated in Figure 6. kach property is plotted in terms of percent 


of the final value. 


c4+u 
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Variation of Properties within the Detonation Zone. 
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At the shock front the pressure is twice the value attained at the end 
of the reaction zone, since the substitution of nQ into Equation 28 


and simplifications yields, to close approximation, the relation 
Pr = (1 + Jv i-n )P} 


The mass velocity u, is always the same fraction of its final 
value as is the pressure P,; therefore, the same curve represents the 
variation of both parameters through the reaction zone. 

The ratio of (ctu) to D exceeds 2 at the shock front and reaches 
unity at the end of the reaction zone, thus it is seen that within the 
reaction zone a rarefaction wave of velocity c+#u will continually over- 
take and decrease the velocity of a discontinuity progressing at a lower 
velocity. At the end of the reaction zone D = ctu, (the Chapman-Jouguet 
condition) and a stable detonation wave is possible. 

The temperature and entropy show sharp increases at the shock 
front and a gradual rise throughout the reaction zone. 

Equation of State Based on Energy Considerations 

In order to improve the agreement between calculated and 
observed values of the detonation velocity, equations of state more 
representative of the extreme conditions confronted have been suggested. 
Jones (27) has developed an equation suitable for use at high pressures 
based upon Bridgman's (28) results with nitrogen at pressures up to 


15,000 atmospheres. The equation takes the form 


V/N 


P= RT £(#,7) - Tore 


where Eo(V/N) is assumed to be the potential energy of interaction of 
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molecules in the gaseous product and N is the number of moles. fhe 


internal energy E, is expressed in the form 
EQ(v/N) = aew “W/N 4 BIV/N) + 


in which the first term represents the repulsive forces and the second 
the attractive forces. The constant C is assigned a value to make the 
minimum energy equal to the heat of vaporization. The function f is 
evaluated from compressibility data and the constants A, B and @ from 
the energy of the normal modes of vibration of molecules in the fully 
excited solid state. 

At lower temperatures, Jones assumes an equation of the virial 
form, the constants being adjusted so that the pressure computed at an 
intermediate value of V by the both relations is the same. Extending 
the interaction effects determined by Bridgeman to higher temperatures 
and densities the total energy of the reaction products can be expressed 


as 


E> S([NgE + ( ZNy-Ng)(Eo+3/2a7)| 


where N; is the number of moles of each molecular species, both gas and 
solid, E; is the energy per mole exclusive of the interaction of 
molecules in the gas phase and Ne is the number of moles of solid 
products. ‘The last term, therefore, accounts for the energy of molecular 
interaction plus the vibration energy of 3/2RT per mole. The energies 

E, can be evaluated from specific heat data at ordinary pressures since 
the interaction effects at high pressures are represented by E,. To 
determine the composition of the reaction products as a function of 


temperature, the activities of the various products are introduced. 
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The activity of each species can be expressed in terms of the composition 
N; and the equilibrium constants for the assumed reactions as a function 
of temperature at normal pressures to which ideal gas conditions are 
applicable. The increase in internal energy of the products must be the 
sum of the chemical energy Q released by the reaction and the work dme 
by the pressure at the detonation front. Thus a knowledge of E, Q and 
the equation of state provides sufficient information to solve the 
detonation front equations. 

Although the above procedure involves several major assumptions 
to make it feasible, the calculation is significant because the results 
tabulated in Table III give an idea of the success obtainable using an 
equation of state fitted to data other than that obtained by detonation 
measurements. 

The Halford-wilson-Kistiakowsky Equation of State 

Information on gas states at the pressures and temperatures 
encountered in detonations can be obtained by the use of the hydro- 
dynamic theory if the variation of detonation velocity is known as a 
function of the initial specific volume (loading density). Since the 
detonation velocities of solid explosives show a definite increase with 
loading density, such measurements provide a suitable means for finding 
the deviations of explosion gases from ideal gas behavior. The detonation 
velocity-loading density relationship is of little value in suggesting 
the correct form of an equation of state, but, once a form is established, 
the relationship is sufficient on the quantitative evaluation of its 
empirical parameters. Unfortunately, the pressure-volume relationship 


is rather insensitive to the form chosen which suggests the measurement 
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of detonation temperature for determining the appropriateness of the 
choice. Since detonation temperatures have not as yet been determined 
to a degree of accuracy sufficient to provide a check, we are forced to 
accept the form which fits the data best in the largest number of cases, 
An equation of state found to be very successful in the cal- 
culation of detonation velocities is one developed by R. S. Halford in 
cooperation with Wilson and Kistiakowsky (17). The equation of state 
has the following form when all the products of reaction are gaseous 
PVM = nRT(14xeP4) (32) 
with KX = Tnyky/T™ vm 


V is the specific volume, n is the number of moles of gas per M grams; 
P, T and R are the pressure, temperature and the gas constant per gram 
mole, respectively; and n; and k; are the number of moles and an empirical 
‘covolume" parameter of the i-th species, respectively; the summation 
extends over all gas species present. Satisfactory results have been 
obtained with the constants “= 0.25 and @ = 0.30. 

In cases in which solid products are formed in the detonation, 
the equation is slightly modified. It is assumed that the volumes of 
gas and solid are aduitive and that the gas obeys Equation 32 independent 
of the presence of solid components. Since thermal expansion and baric 
compression of solids have opposing effects in respect to the volume, and 
since these effects are not known at the elevated temperatures and 
pressures encountered, it is assumed that the specific volume of a solid 
component is independent of temperature and pressure. Thus the specific 
volume of solid components may be assigned values determined under normal 


conditions or may be treated as adjustable parameters. With these 
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assumptions, the equation of state for a mixture of solid and gaseous 
components becomes 


PM(V-V,) = nRTF(x) (33) 


where Vg is the volume of solid constituents per gram of mixture, 


obtained by 


Vs =™a Va * "py Vy + eee (33a) 


Ma? peers are the weight fractions and V,, Vp,... are the specific 


volumes of solid constituents a, b,... F(x) is a function defined by 


F(x) = 1 + xe &* (33d) 
with 

x = k/T™ M(V-V,) (33¢) 

k = Ink, (33d) 


The value of « and @ and all other symbols are the same as those used 
with Equation 32. A theoretical evaluation of the covolume factors 
(c.f. Lquation 33d) is impractical in view of the extremely high density 
encountered in the reaction zone. These factors are therefore determined 
by comparison of the results of a calculation of detonation velocity 
based on an assumed value with the actual observed velocity. As an 
initial approximation only the most abundant detonation products are 
considered and the total covolume constant k of the mixture is obtained 
for several explosives at a number of loading densities. To evaluate 
the specific covolume factors k;; it is assumed that k is an additive 
function of specific values. ‘The individual values are then adjusted 


by an iterative procedure for optimum conformance with the mixture values. 
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Since the original work was completed, other methods of evaluating these 
constants have been suggested and a wider range of explosives has been 
used as a basis, yielding slightly modified values. It has been found 
that the calculation of product composition is much more sensitive to 
the individual values of ky than is detonation velocity. 

The application of the Halford-Wilson-Kistiakowsky equation 
of state to the calculation of detonation velocities and related 


quantities will be illustrated below. 
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The calculation of the detonation velocity and related 
quantities using the Halford-wilson-Kistiakowsky equation of state takes 
into account the effects of gas imperfections and results in an entirely 
acceptable set of parameters in the vast majority of explosives to which 
it has been applied. 


Calculation of Detonation Velocity 
From pure thermodynamics, we obtain the relation 


av), (atl ylasl, fan), Gu, ~ "la, 9 


Differentiating the equation of state, Equation 33, holding T 


constant gives the result 
dap =- Seal Pe g xe vs) (35) 
From Equation 33b, the coefficient of (V-V,) can be represented 


by F(x) and that of P by y(x) where 


y(x) = F(x) + an (36) 
Thus kquation 35 can be written 

ap) = = Py(x = 1 (ap 

(au V-Vg) F(x + (i). (37) 


Substitution into the basic equation for D, Equation 21, yields the 


expression 


D= Vo P)T y Xx) V,-Vs F x) (38) 
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Introducing the equation of state once again gives 


. V 
Dus . ma ¥ ny RT} y (xy )/M (39) 
S 


1 





Equation 38 may be substituted into Equation 4 to obtain an expression 


in terms of the specific volumes (P, is neglected) 


p, = Paya) Fi Wo-Va) (40) 
which on rearrangement becomes 


lL = 1+F(x, ) 
Vy-V, y(x4y Fy 





(41) 


Calculation of the Specific Heat Ratio 


From thermodynamics, one obtains the following equations for 
the specific heat at constant volume and the ratio of the specific 
heats at constant pressure to that at constant volume 


C, = Cy + un{Y-Ys (d2p/aT*) av (42) 
where ca is the ideal heat capacity per M grams and 


yea Bl, /(Gl, © 


Introducing the equation of state in these expressions yields 


Cy = or. + mRX(z(x) = 1) (41) 
and = Y= 1 + mR | 2(x)?*/cyy(x)] (45) 


where z(x) = F(x) -Xx “Ab (46) 
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Calculation of the Temperature 


For an imperfect gas, Equation 16 expresses the application of 


the first law of thermodynamics 
_. + V,-V ; 
M(E] - Ey) = QeCy (1 )-T,)+M : 1 )(aB/aV) 7, dv (47) 
90 


(Ihe reaction is first carried out to the final temperature at a large 
volume, such that the product gases are ideal, and these gases are then 
compressed to their final volume with a resultant energy term due to the 
gas imperfections. ) 

Introducing the equation of state and standard thermodynamic 


relations, one obtains 
M (1 8lax/av) dqdV = n RT, a( F(x )-1) (48) 
oe Ta ghee Aiea 


Equation 41 can be rearranged to the form 
voqv, = FOq) Wy-¥,) (49) 
¥(4) vy 
which when combined with the Hugoniot Equation 6 yields 


E)-E, = BP) F(x, Vy /¥(xq) v7) = dny RT} FX (x, )/y(x,) OM (50) 


Therefore, from Equations 47, 48 and 50, one obtains 


Ty = Cy To-Q/GyenqR& (F( xy )-1)-nyRPA( x) )/2y(x,) ¥ (52) 


solution of the above equations is based on an assumed knowledge of the 
final composition, thus an iterative approach appears desirable. If a 

large number of explosives are to be investigated a two-stage method is 
recommended, This involves the calculation of ideal values of D, T, and 


v based on the assumption that the product gases are ideal and then 
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evaluating the ratio of actual to ideal values using tables prepared 
for this purpose. 
The entropy increase occurring in the detonation zone can be 
estimated if values for the increase in temperature and volume across 
the shock wave are assumed. 
The entropy increase across the shock zone due to an irreversible 


compression to a pressure P and temperature T is given by 


4s = Cy In + | aa (Vo-¥y) oe 


The entropy change due to the irreversible chemical reaction 
can be calculated by proceeding from the initial to the final state by 


hypothetically reversible steps, applying the relation 


T R 
AS = S°pog. - S'explos, #NCpln =t-nRin i (53) 
Prod. Explos pin 1, PS 


in which the subscript 1 refers to the conditions prevailing at the end 


of the reaction zone. 
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COMPOSITION OF THE DETONATION PRODUCTS 


In order to carry out calculations of the detonation velocity 
and other parameters of interest, it is necessary either to calculate 
the equilibrium composition of the reaction products or to make an 
assumption regarding this composition. Before discussing the detailed 
methods used to determine the composition, a few general remarks concern=- 
ing equilibria among possible products appear desirable. 

Dissociation into free atoms and radicals is fortunately not 
of importance at the high pressures encountered in detonation of solid 
explosions. The dissociation constant of H> = 2H at 5000°K has been 
calculated statistically to be 44.7 atm. Since the concentration of 
hydrogen gas rarely exceeds 10% by volume in most explosives, it is found 
that at a pressure 10? atmospheres about 2% of the hydrogen is 
dissociated into atoms. This corresponds to an absorption of heat of 
about 3 kilocalories per kilogram of explosive, whose total heat of 
explosion is of the order of 1000 kcal. per kg. Thus, the final 
temperature is lowered to the extent of 0.3% by the inclusion of this 
dissociation but at the same time, the number of moles of products and 
the pressure is increased. The two effects oppose each other as far as 
detonation velocity is concerned and it is estimated that the net effect 
results in a 0.10% decrease in the velocity. The dissociations of Oo 
and N5 occur to a much lesser extent and can also be ignored in all 
calculations. Although the dissociations of water: H20 = 3 Hp + OH 
and Ho0 = 3 O02 + Ho take place to only a very small extent these 
equilibria are considered in precise calculations. 


Various polyatomic molecules have been detected in gas samples 


\ 
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analyzed after closed bomb detonations. Schmidt (29) has reported small 
but variable amounts of CH), NH CoHo, HCN and CoNo from almost all 
explosives. It is suspected that the majority of these molecules are 
formed after the detonation is complete and while the products are 
cooling. Allowance is made in precise computations for the formation of 
methane and ammonia, but the remaining three are neglected in the 
interest of mathematical simplification. If it is assumed that the 
analytical data represents the true conditions existing in the detonation 
wave, the formation of all these complex molecules has a small but not 
negligible effect on the results of calculation. For example, in TNT 

. where the oxygen deficiency is large and the polyatomic molecules listed 
above are particularly abundant neglecting their formation results in an 
overestimation of the detonation velocity by about 3.5%. 

The formation of free hydrogen in absence of solid carbon is 
considered in precise computations even though the equilibrium constant 
of the "water-gas" reaction, H90 + CO = COotHo ranges from 0.2 at 3000° 
to 0.07 at 5000°K. Since the number of moles of products does not change 
the only effect on the calculations is that due to the evolution of heat. 
The neglect of this equilibria in the absence of solid carbon can result 
in an underestimation of the rate by only a few percent. If, on the 
other hand, free carbon is present, the reaction C+H20 = CO+H2 proceeds 
almost quantitatively. The equilibrium constant of this reaction is of 
the order of 10? atm. at 3000° and 10° atm. at 5000°K. Since the 
fugacity of CO in the reaction zone is usually less than 10° atm., the 
ratio fy,,/ fHD0 is greater than unity as long as carbon is present. 


The assumption that the reaction goes completely to the right results in 
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an overestimation of the detonation rate since the heat absorption is 
overcompensated for by an increase in the number of moles of gaseous 
products. 

As a first approximation results may be obtained using the 
following arbitrary decomposition equations for explosives of the 
general formula, CGHOSN « 

Case A: (q+r/2€S8) 


CHLO.N, = r/2 H20 + (s - q - r/2)CO, + (2q- 8 + r/2)CO + t/2 No 


Case B: (q + r/22 S) 


C HON, = r/2 H90 + (s - r/2)CO (q-s + r/2)C + t/2 No 


These equations express the assumption that all the hydrogen 
is converted to water, that the remaining oxygen then converts carbon 
to carbon monoxide and that additional oxygen, if any, converts carbon 
monoxide to carbon dioxide. The results of a calculation using the 
arbitrary product composition given by the above equations are shown in 
Table III. 

If one assumes that the detonation products are in chemical 
equilibrium at the final temperature pressure and density, then the 
equilibrium relations together with the stoichiometric restrictions 
provide the necessary relations for the calculation of the product 
composition. The thermodynamic equilibrium constants for the significant 
equilibria can be computed from the function - (FO-E3)/T for each of 
the components, together with AES for the equilibria of interest. 
Tables of these quantities for Hoy CO, COz, Ho0, No, OH, Og, and NO 


and for the following equilibria over the temperature interval, 2000° - 
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5000°K. have been published (30). 
Ho + COo = CO + Ho0 
200 = COo + C 
HO = 3Ho + OH 
2H20 = 2H2 + Oo 
Ho0 + aN5 = NO + Ho 
The problem is conveniently approached by the consideration of 
the work content A of M grams of a mixture of gases at a large volume 
v° and temperature T}. From the thermodynamic relations, we have 


A(V®, T]) = Ing wg? + RT, In7t(p$) 2 ~ nkT, (54) 


where nj is the number of moles of the i-th constituent and M; and 
pi are its chemical potential and partial pressure, respectively, in 
the standard state (pure i-th component at T, and unit pressure). 
nis the total number of moles. 


Equation 54 can be written 


A(V°, T}) = Inj ag + RT) [n ln P/n + 2ny in ny — n] (55) 


The work content at a volume V, is obtained by 


A(v°, T)) -M f dpav (56) 


To evaluate the integral requires the introduction of an equation of 

State. Since calculations reported in a later section have been based 
on the Halford-Wilson-Kistiakowsky equation, Equation 33, this will be 
used, To simplify the resulting expression it is convenient to use the 


function F(x) as defined in Equation 33b,. Performing the indicated 
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integration and simplifying one obtains 
A(V,T]) = Bny~P+RT, {ny Inn, + nln [ P,/nF(x,)] 
n(e®* ~ 1) /@ -n} (57) 


The chemical potential is given by 


AL = (2A/ 2n4)0 yin (58) 


and therefore 


Ay = Ay + RT, [ In n; + In(P}/n) + 2.303 hi (xy) + 


2.303 n ky F2(x)/«] (59) 
where [3 (x) ® 0.4343 (e a 1)/@ - log F(x) (60) 
[> (x) = 0.4343 x e OX (61) 


These functions have been tabulated with x as the argument (31). The 
chemical reaction may be represented by 


EWss x = 0 (62) 


where the mole number Ws for the i-th constituent in the j-th reaction 
is taken positive for products and negative for reactants. Thus 
Equation 59 may be written 
Oo 
2.303 F(x) + 2.303 n DWyy ky MZ(x)/k] (63) 


Now, since at equilibrium 2. W54 A; = O and by definition 2 W535 My = 
- RT ln RD we may write 
log MP4 = 2 Wiss log n, - Log G, (64) 


where Ky j is the thermodynamic equilibrium constant for the j-th reaction 
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and the quantity G, is defined by the equation 


J 


~log G, = [20g(P,/n) +3 (x) | 2Ws5 +n [> (x ) 
2 W54 k,/k (65) 


Thus one can now write for the equilibrium constant in terms of 
moles (K,) 


W 
= 77 Ja 
p, 05 = (my) (66) 


Kk, =k 
If there are 1 constituents in the mixture, composed of m different 
atomic species, there will be (1 - m) relations of the form of 
Equation 66 and m stoichiometric relations. A convenient iterative 
procedure and the tabulated functions necessary for the calculation of 
the equilibrium composition of the reaction products are available in 
reference (31). The equilibria usually considered include the five 
equations given above as well as 

3H, + CO = CH, + HAO 


and 3/2 Ho + 4 No = NHs 


A comparison of the final composition of some typical explosives 
based upon the arbitrary decomposition equation and the equilibrium 
considerations is given in Table II below. The composition is given 


as the number of moles of each constituent per mole of explosive. 
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Table II 


INT Py = 1.0 Tetryl 0, = 1.0 Cyclonite Q, = 1.3 


Arbit. | Equil, | Arbit. | Equil. | Arbit.| Equi 





Comparison of product composition, based upon arbitrary and equilibrium 
relations. 
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RESULTS OF THE CALCULATIONS 


Table III is an example of the results obtained by Brinkley 
and Wilson (31) using the method of approach discussed in the previous 
section. The initial and final densities, @, and (1, are in grams/cubic 
centimeter. The detonation velocity, D, and the mass velocity, u,, are 
in meters per second. The temperature is in degrees Kelvin and the 
pressure in kilobars (107 dynes/cm@). The composition is in gram moles 
of the individual constituent per gram mole of explosive. 

A tabulation of the results obtained by Jones (27) using INT 
(°C, = 1.5 is included for comparison. Jones' final composition gives 
larger values for CQp and CH), at the expense of Ho0 and CO. The avail- 
able experimental data on detonation product composition indicate much 
smaller quantities of methane, however, it is entirely possible that 
formation may have occurred during the cooling process, The equilibrium 
of reactions involving CO and COs are also sensitive to temperature changes. 
On the whole, it is reasonable to suppose that the compositions given 
by the Brinkley-Wilson procedure are a closer approach to the actual 
conditions. This is also suggested by the better agreement with the 
observed (32) detonation velocity of 6620 meter/sec. for an initial 
density of 1.5. 

the right-hand column contains results based on the arbitrary 
decomposition equations. It appears that the detonation velocity, which 
is the only property that can be measured accurately, is comparatively 
insensitive to the composition. Although the three methods give results 
of the same order of magnitude, it is felt that the equilibrium approach 


of Brinkley and Wilson is preferable. 
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A comparison of observed and calculated detonation velocities 
of some common explosives is shown in Figures 7, 8 and 9. The observed 
velocities are those reported by the kxplosives Research Laboratory (32), 
Bruceton, Pa. The calculated values are based on the method of Brinkley 
and Wilson with consideration of the seven primary equilibria. The 
overall agreement is considered to be very good. 

Deviations may be due to errors in the observed quantities, 
however, these are usually accurate to within 3 - 4%, or to errors in 
the heat of formation of the explosive or to inadequacies of the method 
due to lack of accurate data in regard to covolume factors and equilibrium 
constants under the extreme conditions. 

The heats of formation of the explosives considered in Figures 


7, 8 and 9 are based on the values reported by Schmidt (33) in 1924. 








Table IIT 


Calculated Values of the Detonation Velocity 


and of the Properties of the Detonation Wave 
TINT 

% tt 
@® Wo li 22 23) 1k See oe bee 
CP, 1.347 1.464 1.578 1.690 1.801 1.911 2.019 1.94 
D 4870 5190 5540 5920 6340 6820 7290 7720 7210 
Uy 1250 1290 1330 1378 1480 eso 510 W750 
Ty 3260 3270 3270 3240 3210 3190 3170 3460 2880 
Py 61.0 73.6 86.6 10590 1250) 15396 G7yes 205 173 
Composition 
ie | 09 05D tC 0 0 0 0 0 
eo, “16 “21 .25 ep 83 6 See 2p Oe 
CO 3.36 3.22 3.08 2.92 2.75 2:58 2542 06 35 
HO 382.260 «62.36 «62.420 «20460 2048 2048 22K 110 25 
No lib lek? 1.48 WAG 1649 Ves ieee 190 1:5 


NO QO QO 0 0 Ol 02 002 Q Q 
CH, .O1 0 0 QO QO 0 Q se O 
NH3 .O? 06 04 02 Ol Ol OL 0 QO 


C 3088 3057 mT eeTy)0=— 3.87) | ae 0m 63 BR 95 





* Based on Jones! equation of state. 
#* Based on arbitrary decomposition equations. 
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APPENDIX A 


RELATION OF DETONATION VELOCITY TO c + u AT VARIOUS POINTS 
ON THE RANKINE-HUGONIOT CURVE 


In a previous section, it was shown that at the point Bb, 
Figure 3, the detonation velocity was equal to c + u in accordance with 
the Chapman-Jouguet condition. It is of interest to determine the 
relation at points other than B. Figure 3 is reproduced below for 


convenient reference. 





Generalized Rankine-Hugoniot Curve 


It is well established that the sonic velocity is related to 


the slope of an isentropic PV curve by 


Re22 | 
Cc, = V -—— 
a : | OVi/s 


From Equations 4 and 5, we have the relation 


2 2 Poetic 
il 1 Vo-Vy 


PoVo being evaluated at the point A and PV) on the R-H curve. From 


these two equations, it appears that the difference between of and 
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(p=u, )* at various points on the R-H curve is 


of - (D-u, )* = vy (- en - eae (1) 
ov} 8 VorVo 


The slopes of the isentrope and R-H curves at any point (P, V) are 


related by the equation 
(25) = (22) (24) E., 
OVip OV], IS }y OVipl R 


wnere the subscript R indicates differentiation along the R-H curve. 


Rearranging and introducing the Maxwell relation cA ae 2 
we have 


- 22) -(22) =. (22) (28 2 
IVs OVip (34), WV] R 2) 
From the energy Equation 6 for points on the R=H curve, we have 
sou (Bh 
Ay. = Vi-V “Ay - P, +P 
SH) 2 Worry) (Buy) 3 rte) 


Substituting for =| its equivalent given by the first and second 
R 





eV} 
laws of thermodynamics yields 
( =) ‘ Vo-¥,. | Pit, oP} 
Min 2% «| VOW, VOM, (3) 


The sign of quantity on the left-hand side of this equation depends upon 
the relative slopes of the line joining the initial and final states and 
the tangent to the curve at the point representing the final state. Two 
extremes are shown in the figure below. Note that at point B the two 


lines coincide, 











Figure 10 


Slope Relations on the R=-H Curves. 





ave 
negative above B and positive below B. It is certainly reasonable to 


2 
If the R-H curve is concave vpwara (2 Bho, the difference in slopes is 


expect that the R-H curve should have this type curvature for the 
majority of normal fluids. Solving Equation 3 for the quantity within 
the square brackets and subtracting from Equation 2, evaluated at the 
point (P)V,) permits the velocity difference written in Equation 1 to 
be expressed in terms of variables of known sign. That is 
ont -(B[ A (ale 
The terms within the square brackets are always positive as is Ve; 
therefore the right-hand side is positive for points on the R=H curve 
above B (where ( 95,/ 9V})p is negative) and negative for points below 
B, that is where ( 9S,/ OV) )p is positive. Hence, as was to be shown, 
D is less than c)+u, when the final state (P,V,) is above B and is greater 


than 6) Pu) when the final state is below B, 
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APPENDIX B 


SECONDARY EFFECTS 


The theoretical treatment of the detonation of solid explosives 
was most conveniently accomplished by neglecting certain effects such as 
charge diameter and unstable detonation velocities that must be considered 
in any practical application of the theory. Fundamental investigations 
of these effects and their ramifications have been made and several 
theories have been proposed in an effort to explain the observed 
phenomena. Since a detailed description of the proposed theories sheds 
little light on the fundamentals of the detonation process, it is believed 
advisable, for the purpose at hand, to regard them as "second-order" 
effects and to devote to such effects a very limited treatment. 

The effect of charge diameter on the detonation velocity has 
been observed since the earliest experiments were conducted. This effect 
is observed as a decrease in the detonation velocity when an explosive is 
confined in a tube of small diameter. Fortunately, the variation ceases 
wnen a certain reasonable diameter is reached, thus the investigation of 
other parameters can be carried out without the introduction of an 
additional variable. By the same token, few investigators have under- 
taken detailed studies of its cause and until recently a generally 
accepted explanation was not available. A theory proposed by Jones (34) 
accounts for the observed variation on the basis of lateral losses. 
During the finite reaction time, the high pressure with the reaction zone 
is assumed to cause an appreciable lateral expansion with a resultant 


drop in pressure, temperature and detonation velocity. In other words, 
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not all of the energy of the explosive is available to maintain the ideal 
detonation velocity due to lateral dissipation. It can be shown by a 
detailed mathematical analysis that the effect becomes significant when 
the radius of the explosive charge is of the same order of magnitude as 
the reaction zone. 

Another phenomenon observed frequently in the detonation of 
solid explosives is the transition from an unstable detonation velocity 
(one that is either too high or too low) to the stable value. Whenever 
an explosive charge is initiated by a priming charge having a different 
detonation rate, this transition must occur. A theory of the building-up 
rate as a function of the reaction zone length has been proposed (7) and 
leads to the conclusion that the ideal velocity is approached asymtotically, 
half the difference between the initial and final value being covered 
in a few reaction zone lengths. 

If a detonation wave experiences such severe losses that the 
wave can no longer propagate at the ideal velocity, the phenomenon 
referred to as "fading" is observed. As in the case of transition from 
one velocity to another, the length of the reaction zone is a critical 
feature and theoretical treatments of the phenomenon are based upon the 
cumulative effects resulting from increased reaction zone lengths. 

In summary, it might be observed that if one possessed sufficient 
information concerning the chemical reaction kinetics to predict the 
length of the reaction zone under any conditions of pressure, temperature, 
etc., a more rational approach to many of the technical problems con- 


fronting the explosives engineer could be made. 
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